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Building on a fundamental insight due to Avramidi, we provide a system of transport equations 
for determining key fundamental bi-tensors, including derivatives of tiie world-function, a{x,x'), 
the square root of the Van Vleck determinant, A^''^ {x,x'), and the tail-term, V{x,x'), appearing 
in the Hadamard form of the Green function. These bi-tensors are central to a broad range of 
problems from radiation reaction to quantum field theory in curved spacetime and quantum gravity. 
Their transport equations may be used either in a semi-recursive approach to determining their 
covariant Taylor series expansions or as the basis of numerical calculations. To illustrate the power 
of the semi-recursive covariant series approach, we present an implementation in Mathematica which 
computes very high order covariant series expansions of these objects. Using this, a moderate 
laptop can, for example, calculate the coincidence limit [07(2;, a;)] and V{x,x') to order {a°')^° in a 
matter of minutes. Results may be output in either a compact notation or in xTensor form. We 
also numerically integrate the transport equations along null geodesies in Nariai and Schwarzschild 
spacetimes. 



I. INTRODUCTION 

In a recent paper [T] we presented methods for obtaining coordinate expansions for the (tail part of the) retarded 
Green function in spherically symmetric spacetimes. By using computer algebra to obtain high order Taylor series 
(of order (Aa;")^'^), and applying the theory of Pade approximants we were able to obtain accurate expressions in 
remarkably large regions. Using these expressions, we were able to present the first complete matched expansion 
calculation of the self-force in a model 'black hole' spacetime, the Nariai spacetime [2 , and are currently applying 
the method to Schwarzschild spacetime. Our ultimate goal in this programme is to work in more general spacetimes, 
especially Kerr spacetime. A key component of the matched expansion approach is knowledge of the Green function 
for points close together (i.e. in a quasilocal region). As we move away from specific symmetry conditions, we can no 
longer rely on methods based on a special choice of coordinates in the construction of our quasilocal solution and are 
led instead to consider other techniques such as transport equations and covariant expansion methods. 

Covariant methods for calculating the Green function of the wave operator and the corresponding heat kernel, briefly 
reviewed in Sec. below, are central to a broad range of problems from radiation reaction to quantum field theory in 
curved spacetime and quantum gravity. There is an extremely extensive literature on this topic; here we provide only 
a very brief overview referring the reader to the reviews by Vassilevich 3 and Poisson ^ and references therein for a 
more complete discussion. These methods have evolved from pioneering work by Hadamard [5j on the classical theory 
and DeWitt [BJ |7] on the quantum theory. The central objects in the Hadamard and DeWitt covariant expansions 
are geometrical bi-tensor coefficients a^^ {x,x') which are commonly called DeWitt or HaMiDeW coefficients in 
the physics literature. These coefficients are closely related to the short proper-time asymptotic expansion of the 
heat kernel of an elliptic operator in a Riemannian space and so are commonly called heat kernel coefficients in the 
mathematics literature. Traditionally most attention has focused on the diagonal value of the heat kernel K^a{x, x; s), 
since the coincidence limits a^A{x,x) play a central role in the classical theory of spectral invariants [8] and in the 
quantum theory of the effective action and trace anomalies [S]. By contrast, for the quasilocal part of the matched 
expansion approach to radiation reaction |10l 111] we seek expansions valid for x and x' as far apart as geometrical 
methods permit. 

The classical approach to the calculation of these coefficients in the physics literature was to use a recursive 
approach developed by DeWitt [J in the 1960s. Although these recursive methods work well for the first few terms in 
the expansion [T^[T^, and may be implemented in a tensor software package [TJ], the amount of calculation required 
to compute subsequent terms quickly becomes prohibitively long, even when implemented as a computer program. 
An alternative approach, more common in the mathematics literature, is to use pseudo-differential operators and 
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invariance theory P] , where a basis of curvature invariants of the appropriate structure is constructed [TS] and then 
their coefficients determined by expHcit evaluation in simple spacetimes. However, here too, the size of the basis grows 
rapidly and there seems little prospect of reaching orders comparable to those we obtained in the highly symmetric 
configurations previously studied. 

An extremely elegant, non-recursive approach to the calculation of HaMiDeW coefficients has been given by 
Avramidi |16[ 117] . but as his motivation was to study the effective action in quantum gravity he was primarily 
interested in the coincidence limit of the HaMiDeW coefficients, while in the self-force problem, as noted above, we 
require point-separated expressions. In addition, Avramidi introduced his method in the language of quantum me- 
chanics; quite distinct from the language of transport equations, such as the Raychaudhuri equation, more familiar to 
discussions of geodesies among relativists. In this paper we present Avramidi 's approach in the language of transport 
equations and show that it is ideal for numerical and symbolic computation. In so doing we are building on the work of 
Decanini and Folacci |18j who wrote many of the equations we present (we indicate below where we deviate from their 
approach) and implemented them explicitly by hand. However, calculations by hand are long and inevitably prone to 
error, particularly for higher spin and for higher order terms in the series and are quite impractical for the very high 
order expansions we would like for radiation reaction calculations. Instead, we use the transport equations as the 
basis for Mathematica code for algebraic calculations and C code for numerical calculations. Rather than presenting 
our results in excessively long equations (our non-canonical expression for 07(2:, x) for a scalar field contains 2 069 538 
terms!), we have provided these codes online [l9l [20]. 

In Sec.|TT] we provide a brief review of Green functions, bi-tensors and covariant expansions, outlining the relations 
between the classical and quantum theories. 

In Sec. HI we detail the principles that we consider to encapsulate the key insights of the Avramidi approach and 
use these to write down a set of transport equations for the key bi-tensors of the theory. These provide an adaptation 
of the Avramidi approach which is ideally suited to implementation on a computer either numerically or symbolically. 

In Sec. |IV| we describe a semi-recursive approach to solving for covariant expansion and briefly describe our 
Mathematica implementation of it and its interface to the tensor software package xTensor. 

In Sec. |V] we present a numerical implementation of the the transport equation approach to the calculation of 
V{x,x') along a null geodesic. 

Given our motivation in studying the radiation reaction problem we shall phrase all the discussion of this paper 
in 4-dimensional spacetime. The reader is referred to Decanini and Folacci [ISJ for a discussion of the corresponding 
situation in spacetimes of more general (integer) dimension. We do note however that the HaMiDeW coefficients are 
purely geometric bi-tensors, formally independent of the spacetime dimension. 

Throughout this paper, we use units in which G = c = 1 and adopt the sign conventions of [3T]. We denote 
symmetrization of indices using brackets (e.g. {aP)) and exclude indices from symmetrization by surrounding them 
by vertical bars (e.g. (a|/3|7)). Roman letters are used for free indices and Greek letters for indices summed over all 
spacetime dimensions. 



II. A BRIEF REVIEW OF GREEN FUNCTIONS, BITENSORS AND COVARIANT EXPANSIONS 



A. Classical Green functions 



We take an arbitrary field Lp"^{x) where ^ denotes the spinorial/tensorial index appropriate to the field, and consider 
wave operators which are second order partial differential operators of the form |17j 

V^B = 5^Bia-m^) + P^B (2.1) 

where □ = ^"''VqV^ and Vq is the covariant derivative defined by a connection A'^Ba'- '^a^^ = da'-P^ + A'^Baf^ , 
m is the mass of the field and P^b{x) is a possible potential term. 

In the classical theory of wave propagation in curved spacetime, a fundamental object is the retarded Green function, 
Grct^c {x^x'). It is a solution of the inhomogcncous wave equation, 

2?^sG,et^c' {x. x') = -4^<5^c"5 (2:, x') , (2.2) 

with support within the past light-cone of the field point. (The factor of An is a matter of convention, our choice here 
is consistent with Ref. ^.) Finding the retarded Green function globally can be extremely hard, however, provided x 
and x' are sufficiently close (within a causal domain) , we can use the Hadamard form for the retarded Green function 
solution [5, .22J , which in 4 spacetime dimensions takes the form 

G,et^B' (x, x') = e. (x, x') {U^B' (x, x') 5 {a {x, x')) - V^B' (x, x') 9 {-a (x, x'))} , (2.3) 
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where 0- (a;, x') is analogous to the Heaviside step- function, being 1 when x' is in the causal past of x, and 
otherwise, 5{x,x') is the covariant form of the Dirac delta function, U^^ {x,x') and V^^ {x,x') are symmetric 
bi-spinors/tensors and are regular for x' — > x. The bi-scalar a{x,x') is the Synge [1] world function, which is e qual 



to one half of the squared geodesic distance between x and x'. The first term, involving U b' {x,x'), in Eq. (2.3 1 
represents the direct part of the Green function while the second term, involving V^b' {x,x'), is known as the tail 
part of the Green function. This tail term represents back-scattering off the spacetime geometry and is, for example, 
responsible for the quasilocal contribution to the self-force. 

Within the Hadamard approach, the symmetric bi-scalar V"^^ {x, x') is expressed in terms of a formal expansion 
in increasing powers of a |18j : 

oo 

V^^' {x, x') = F/^' {x, x') {x, x') (2.4) 

The coefficients U^^ and Vr"^^ are determined by imposing the wave equation, using the identity a-aO"" — 2a — 
a-a'C'"' , and setting the coefficient of each manifest power of equal to zero. Since b' is symmetric for self-adjoint 
wave operators we are free to apply the wave equation either at x or at x'; here we choose to apply it at x' . We find 
that U"^^' (x, x') = A^/^ {x, x') g^^' {x, x'), where A (x, x') is the Van Vleck-Morette determinant defined as |5 

A {x, x') = - [-9 {x)]-^'^ det {-a,^p, (x, x')) [-g (x')]"'^' - det (x, x') a'",, (x, x')) (2.5) 

with 5" a {xtx') being the bi- vector of parallel transport (defined fully below). In making this identification we have 
used the transport equation for the Van Vleck-Morette determinant: 

cr'"ValnA = (4- Dcr). (2.6) 
The coefficients V^^ {x,x') satisfy the recursion relations 

^;a'(^-l/2^AB')^, + (r + 1) A-l/2yAB' ^ L^-m-J^B' ^^yAC^ ^ g (2.7a) 

2r 

for r e N along with the 'initial condition' 

a^"'(A-l/2^AB')^^,^^-l/2^AB'^ l^-l/2pB'^,(^l/2^AC') ^ q (2.7b) 

These are transport equations which may be solved in principle within a causal domain by direct integration along the 
geodesic from x to x'. The complication is that the calculation of V^^ requires the calculation of second derivatives 
of VJ,-? the geodesic; we address this issue below. 

Finally we note that the Hadamard expansion is an ansatz not a Taylor series. For example, in deSitter spacetime 
for a conformally invariant scalar theory all the K-'s are non-zero while V = Q. 

B. The quantum theory 

In curved spacetime a fundamental object of interest is the Feynman Green function defined for a quantum field 
ip'^{x) in the state l^*) by 



Gt^ (x,x') = iT (^|(^^(x)(^^ (x')l^) 

where T denotes time-ordering. The Feynman Green function may be related to the advanced and retarded Green 
functions of the classical theory by the covariant commutation relations [7j 

Gf^^'(x,x') = ^ (Gfdf (x,x') + G^/(x,x')) + '-m\x)^^\x') + ^''\x')0^{xm. 

The anticommutator function {^\(p^{x)Lp^ (x') -I- (p^ {x')(p^{x)\^) clearly satisfies the homogeneous wave equation 
so that the Feynman Green function satisfies the equation 

V^BGf^'c' {x,x')^-5^c'5{x,x'). 
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Using the proper-time formalism fT", the identity 

oo 

i / ds e"*^" exp(isa;) = , (e > 0), 

J x + ie 



allows the causal properties of the Feynman function to be encapsulated in the formal expression 

oo 

G{'^C'{x,x')^i [ dse""exp(isI?)^B(5^c"5(x,x') 



where the limit e — > 0+ is understood. The integrand 

K-^C' {x, x'; s) = exp{isV)^BS^c'S{x, x') (2., 
clearly satisfies the Schrodinger/heat equation 

1 dK^c 



i ds 



-{x,x'-s)^V''bK''c'{x,x'-s) (2.9) 



together with the initial condition b'{xtx' ;()) — 5^ b'{x,x'). The trivial way in which the mass m enters these 
equations allows it to be eliminated through the prescription 

K^C'{x,x'-s) = e-''^"'Ko^c'{x.x'-s), (2.10) 
with the massless heat kernel satisfying the equation 

1 dKfi^ ■ 



ds 



-{x,x';s) = {5''Bn + P''B)KQ"c'{x,x'-s) (2.11) 



together with the 'initial condition' Kq'^ B'ix^x' = 5^ b'5{x,x'). 

In 4-dimensional Minkowski spacetime without potential, the massless heat kernel is readily obtained as 

B'{x,x' \s) = - — ^— -expf — —^S^B' (flat spacetime). (2-12) 
(47rs)^ V 2is/ 

This motivates the ansatz [7J that in general the massless heat kernel allows the representation 

K^^B'{.x,x'-s) ^ j^exv[-^^ ^"\x,x')n^B'{x.x'-s) , (2.13) 
where ri^5/(x,x'; s) possesses the following asymptotic expansion as s ^ 0+: 

oo 

n^B'{x,x';s) ^Y.''rB'{x,x'){tsy , (2.14) 



r=0 



with aQ^B'{x,x) — S'^B' and b'{x, x') has dimension (length)^^*". The inclusion of the explicit factor of A^/^ is 
simply a matter of convention; by including it we are following DeWitt, but many authors, including Decanini and 
Folacci, choose instead to include it in the series coefficients 

Ar'^B'{x,x') = Ai/^a/s'(a;,x'). (2.15) 

It is clearly trivial to convert between the two conventions and, in any case, the coincidence limits agree. 

Now, requiring our expansion to satisfy Eq. (2.11 1 and using the symmetry of Q^b'{x,x'; s) to allow operators to 
act at x' , we find that H,^ b'{x, x' ; s) must satisfy 



1 on 



AB' 



i ds 



Inserting the expansion Eq. (2.14), the coefficients a^^ {x,x') satisfy the recursion relations 

a^-'a^f + (r + 1) a^f - A-'/'{6^' c'D + P^' c) (A^/^a/^') = (2.16a) 
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for n G N along with the 'initial condition' 



■,a' AB' n 



(2.16b) 



with the implicit requirement that they be regular 

Comparing (2.7 1 and (2.161, one can see that the Hadamard and (mass-independent) HaMiDeW coefficients are 
related for a theory of mass m by 



with inverse 



-,2\r-fe+l 



k=0 ^ ' 



ak^B'{x,x') 



In particular, 



^+^^^(mY-'=y,V(a:,a:')+ ^™ ^ 



2\r+l 



k=0 



{r-ky. 



(r + 1)! 



(-,-') = ^^;^y?^(-l)^+^a.+/B'(x,x'). 



2'-+ir! 

These relations enable us to relate the 'tail term' of the massive theory to that of massless theory by 



(— 2m^(Tr/^ 



where Jr{x) are Bessel functions of the first kind. 



C. Classical Approach to Covariant Expansion Calculations 

The Synge world-function, a{x,x') is a bi-scalar, i.e. a scalar at x and at x' defined to be equal to half the square 
of the geodesic distance between the two points. The world-function may be defined by the fundamental identity 

cr„cr" 2(7 = CTc-'O-"', (2.17) 

together with the boundary condition lim a(x,x') — and lim aabix,x') = gab{x). Here we indicate derivatives at 
the (un-)primed point by (un-)primed indices: 

CT^'eeW cra = VaCr CT°' EE V'ct CTa' = Va'CT. (2.18) 

(T° is a vector at x of length equal to the geodesic distance between x and a;', tangent to the geodesic at x and oriented 
in the direction x' — s- a; while a" is a vector at x' of length equal to the geodesic distance between x and x' , tangent 
to the geodesic at x and oriented in the opposite direction. 
The covariant derivatives of a may be written as 

a''{x,x')^{s- s')u'' a'''{x,x') = {s' -s)u''' (2.19) 

where s is an afHne parameter and is tangent to the geodesic. For time-like geodesies, s may be taken as the proper 
time along the geodesic while u"" is the 4-velocity of the particle and 

a{x,x')^-\{s-s')\ (2.20) 

For null geodesies, is null and (j{x, x') — 0. 

Another bi-tensor of frequent interest is the bi- vector of parallel transport, gab' defined by the equation 



ffab'iQO"" = = gab';a'<7' 



(2.21) 
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with initial condition lim gab'{x, x') = gab{x). From the definition of a geodesic it follows that 

x' — *x 

daa'O-"' = -ctq and gaa'cr" = -aa' (2.22) 

Given a bi-tensor at x, the parallel displacement bi- vector allows us write as a bi-tensor at x' , obtained by 
parallel transporting Ta along the geodesic from x to x' and vice-versa, 

T^g\, = T^.g-^'^ = T,. (2.23) 

Any bi-tensor Ta-^^...a^a{---a'^ may be expanded in a local covariant Taylor series about the point x: 

°° ^ (—1)" 

-^ai---a^a^---a^ (-^i ) — ^ ^ j ^ai---aj„a^---a^ ai-'-a^ ('^)'-^ ^ ■ ■ • (J (2.24) 

fc=0 

where the ta;^...a^a'-^^---a'^ Qi -ctfc ^'^^ the coefficients of the series and are local tensors at x. Similarly, we can also expand 
about x': 

Ta\---a„^a\---a'^ (x, X ) ^ ^ j ^ai ■■■a,„a^ ---a^ Q'^---aj.(*^)^^''''^'^ (2.25) 

fe=0 

For many fundamental bi-tensors, one would typically use the DeWitt approach to determine the coefficients in 
these expansions as follows: 

1. Take covariant derivatives of the defining equation for the bi-tensor (the number of derivatives required depends 
on the order of the term to be found). 

2. Replace all known terms with their coincidence limit, x x' . 

3. Sort covariant derivatives, introducing Riemann tensor terms in the process. 

4. Take the coincidence limit x' ^ a; of the result. 

This method allows all coefficients to be determined recursively in terms of lower order coefficients and Riemann 
tensor polynomials. Although this method proves effective for determining the lowest few order terms by hand and 
can be readily implemented in software, it does not scale well and it is not long before the computation time required 
to calculate the next term is prohibitively large. It is therefore desirable to find an alternative approach which is more 
efficient and better suited to implementation in software. In the next section, we will describe one such approach 
which proves to be highly efficient. 



III. AVRAMIDI APPROACH TO COVARIANT EXPANSION CALCULATIONS 



Avramidi |16j has proposed an alternative, non-recursive method to the calculation of the coefficients of covariant 
expansions of fundamental bi-tensors. Translated into our language, this approach emphasizes two fundamental 
principles when doing calculations: 

1. When expanding about x, always try to take derivatives at x' . The result is that derivatives only act on the 
CT^'s and not on the coefficients. 

2. Where possible, whenever taking a covariant derivative, Va', contract the derivative with cr" 

Applying these two principles, Avramidi has derived non-recursive expressions for the coefficients of covariant expan- 
sions of several bi-tensors. As Avramidi's derivations use a rather abstract notation, we will now briefly review his 
technique in a more explicit notation. We will also extend the derivation to include several other bi-tensors and note 



that equations (3.11), (3.13), (3.15), (3.17), (3.34), (3.35), (3.46) and (3.49) were previously written down and used 



by Decanini and Folacci [18 

Throughout this section, we fix the base point x and allow it to be connected to any other point x' by a geodesic. 
In all cases, we expand about the fixed point, x. 
Defining the transport operators D and D' as 



D = a"V, 



(3.1) 



7 



we can rewrite Eq. (2.171 as 

{D-2)<7 = {D'-2)a = 0. (3.2) 

Differentiating the second of these equations at x and at x' , we get 

(£''-1)ct° = (£>' - 1) a"' = 0, (3.3) 

which, defining 
can be rewritten as 

a^ = il\,cj'^' a'''=r'„,a«'. (3.5) 



Finally, we define 7° , , the inverse of 



and also introduce the definition 



r\^{ri\,)-'- (3.6) 



^ <y\. (3.7) 



We will now derive transport equations for each of these newly introduced quantities along with some others which 
will be defined as required. Many of these derivations involve considerable index manipulations and are most easily 
(and accurately) done using a tensor software package such as xTcnsor 23J . 

The transport equations of this section may be derived in a recursive manner, making use of the identities 

D'{(Ta[...a'^a' = V^' _^^iD' 0-a[...a'J " a' a' O a'^...a'^ + a'^a' ^^a' c' . . .a'^ + h Cr" -R'' a'^a' ^^a'<^a\...d 

(3.8) 

^'(^^'ai.-.a'J = V^(i?W,...aJ - 77'a'V"'aai...a.^. (3.9) 

and its generalisation, given below. This method is particularly algorithmic and well suited to implementation on a 
computer, thus allowing for the automated derivation of a transport equation for an arbitrary number of derivatives 
of a bi-tensor. 



A. Transport equation for ^"j,/ 



Taking a primed derivative of the second equation in (3.5 1, we get 



(3.10) 



We now commute the last two covariant derivatives in the first term on the right hand side of this equation and 
rearrange to get: 



(3.11) 



B. Transport equation for rfy 



Taking a primed derivative of the first equation in (3.5), we get 

V b' ~ V a'b'(^ + V h' 



(3.12) 



In this case, since cr" is a scalar at x' , we can commute the two primed covariant derivatives in the first term on 
the right hand side of this equation without introducing a Riemann term. Rearranging, we get: 



(3.13) 
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C. Transport equation for 7° 



Solving Eq. (3.131 for ,, and using (3.6), we get 



Next, substituting Eq. (3.141 into Eq. (3.11) and rearranging, we get a transport equation for 7°, 



(3.14) 
(3.15) 



D. Equation for A^^ 



Differentiating Eq. (2.17) at x and x', we get 



which is easily rearranged to give an equation for A";,: 



(3.16) 
(3.17) 



E. Transport equation for a" ^>^i 



Applying the identity (3.8) to (3.111 and simplifying the resulting expression, we get 



V / b'c' c' a o ' h' ol' c' ol' b'c' a.'b'(3'\c' 

ol'p'o' c' a'p'd h' h' p' c' ex' 



(3.18) 



F. Transport equation for (T°'y^i 



Applying the identity (3.9) to (3.11) and simplifying the resulting expression, we get 



(3.19) 



G. Transport equation for ct" 



Applying the identity (3.8) to (3.18) and simplifying the resulting expression, we get 



V -' b'c'd' a'b'c' a' a'b'd' c' ol' d d' h' a'b' c' d' ol' d b' d' a' d' b'd 

+ ^"'.'^"'.v'.' + R^'oL'p'c'R'^'d'yb''^''''^'' + R^'oL'P'b'R'^'d'.'d^''^''' + R'^'oL'P'd'R'^'c'b'^''^'' 



R\'oL'd'R\','d-''^''^ -R\'.'o'R''\','d'^''<^'^ +R\'b'y-.dd''^''<^'^ +R''\' 



/3' 7' 



a' 7' 



''/3'd:d' 



r/3' a' 



^'^'a'.'c';.'^"'^"'.' - R'''oL'P'b';d''^''''^'''d " R''' ol' 0' b' ■,d d' + R'^' d 0' d' <^'' b' .' + R""' b' P' d' d .' 

+ R'''b','d-''^'''d'oL' - R'''oL'p'd'<^''^'''b'd - R''\','d^^'<^'''b'd' - R''\'[3'b'^''^'''dd' = (3.20) 



H. Transport equation for (y°'y^i^i 



Applying the identity (3.9) to (3.18) and simplifying the resulting expression, we get 



/ T^' 1 \ I _OL I a OL I a „CL I a _ql i _a „Oi i a _a 

{U -1)0- b/c'd'+O^ OL'b'd^^ a'b'd'^ d+^ OL'dd'(^ b'+^ OL'b''^ dd'+^ OL'd^ b'd'+^ Q'd'f^ b'd 

+ 'y\"y'''b'dd' + R'''b'P'd;d'''^'^\' + R'''b'p'd'^''^''d'o.' + R''\,i,.,,aP'a\i^, + aV„- = (3.21) 



I. Transport equation for g^,^ 
The bi-vector of parallel transport is defined by the transport equation 



(3.22) 



J. Transport equation for Qab 



Let 



Aabc = Qb" 9c gaa'-p' 

Applying D' and commuting covariant derivatives, we get a transport equation for Aabc'- 

-I'd (3' 



D'Aabc + Aabai^ ^,90,'^ 9c' + ffa^'ff/ 5 -Ra'/3'7'<5' 0"*' = 



(3.23) 

(3.24) 



Let 



K. Transport equation for gab';c 



Babe = 9b 9af3 



Applying D' and rearranging, we get a transport equation for Bafi-y 

I3'9c 



D'Babc = -Aabo.V"0'9c'^' 



(3.25) 
(3.26) 



Applying D' to gj' -c'd'' 



L. Transport equation for 



D 9a -.c'd' ^'^ 9a -.c'd'c 



Commuting covariant derivatives on the right hand side, this becomes 



(3.27) 



D'9:';c'd' - (Sa"'; + R'\' 0' d' 9 a'^' ,c' + c.' 0' C 9 a^' -^d' " R'^' C 0' d' 9 a' + ^"'o'/J'c^.'S/') ■ (3-28) 

Bringing inside the derivative in the first time on the right hand side, and noting that ■0''^^ ^ 0' this then 
yields a transport equation for gj' : 

D'9a';c'd' = -<^^ c'9a^'-0'd' - '^^ d'9a';0'c' " c'd'9a';P' 

+ R'^'^'d'^^'ga"-'-. + R'\.0'c'^''9a''':d' ' ^"c'.'.'^^'s/' + R'' c.' B' C ■d''^'' 9 a^^' ■ (3-29) 



M. Transport equation for = In A^''^ 

The Van Vleck determinant, A is a bi-scalar defined by 

A (x, x') ^ det [A"'^,j , A"'^, ^ -5"'„aV = -ff^'a^V 



(3.30) 



By Eq. (3.13), we can write the second equation here as: 



(3.31) 
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Since D'q" = „,a^ = 0, we can rewrite this as 



A" — D'A" + A" P 
Introducing the inverse (A"-'^)" ^, and multiplying it by the above, we get 



4 = r„- + ^'(lnA) 



(3.32) 



(3.33) 



where we have used the matrix identity Jlndet M = TrM ^(5M to convert the trace to a determinant. This can also 
be written in terms of A^/^: 



^'(lnAi/2)=l(4-r'„, 



(3.34) 



N. Transport equation for the Van Vleck determinant, A^/'^ 

By the definition of C, the Van Vleck determinant is given by 

Ai/2 = e?, 

and so satisfies the transport equation 

2 V ^ " 



(3.35) 



(3.36) 



O. Equation for C^'^l'^ D{lS}l'^) 
Defining r — A^^/^_D(A^/^), it is immediately clear that 

r = A-i/2i5(Ai/2) = 



(3.37) 



P. Equation for 1\-^I'^D' {IS}!'^) 

Defining r' — A^-^/^I?'(A^/^), it is immediately clear that 

r' = A-i/2i^'(Ai/2) = D'C 



(3.38) 



Q. Equation for V^'A 
To derive an equation for Va'A, we note that 



A = det 

and make use of Jacobi's matrix identity 



det [ryV]det 



9a 



d (det A) = tr (adj (A) dA) 

= (det A) tr (A^MA) 



where the operator d indicates a derivative. Applying (3.401 to (3.39), we get an equation for Va'A 



(3.39) 



(3.40) 



Va'A = -A 



(3.41) 
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R. Equation for D'A 



Applying Jacobi's identity twice, together with d(A ^) — —A ^(dA)A ^, we find an identity for the second 
derivative of a matrix: 



d^ (det A) = (det A) (tr (A^MA) tr (A^MA) - tr (A^MAA-MA) + tr (A^M^A)) 



(3.42) 



Using this identity in Eq. (3.391, we get an equation for D'A, 

- {r'^'^Wn'','^'^/) + {g^'^-g/,/) + (3.43) 



(a 8' 3 a' :u' 

9a' 9^ ,^,,9f}.9ii 



S. Equation for D'A^/^ 



Noting that 



□'Ai/2^ f^A-i/^A.^^V = -A-i/2q'a_ Ia-^/^a^a^'a..,, 
1 2 / 2 4 



it is straightforward to use Eqs. (3.411 and (3.431 to find an equation for D'A^/^ 



□'Ai/2 = -Ai/2 
2 



(3.44) 



- {l^'a-"-, + {9a'-9/,/) + (7":-"„,^;,)] (3.45) 



T. Transport equation for Vb 



As is given in Eq. (2.7b), Vq satisfies the transport equation 



or equivalently 



{D' + 1) + Ivr' {e\, - 4) + ip^'c'(AV2g^c') ^ 0^ 



(D' + 1) (a-i/2^ab'^ ^ iA-i/2l?^'c'(Ai/2g-4c") ^ q 



In particular, for a scalar field: 

{D' + l)Vo+ \vo - 4) + ^(D' - - ei?')Ai/2 = 0. 



(3.46) 



(3.47) 



(3.48) 



U. Transport equations for Vr 



As is given in Eq. (2.7a I, Vr satisfies the transport equation 



or equivalently 



iD' + r + 1) 1//^' + iy/^' [e\, - 4) + ^P^'c'KlT = 0. 



(D' + r+l) fA-i/V/^') + lA-i/2pB'^,^AC' ^ 0^ 
v / 2r 



(3.49) 



(3.50) 



Comparing with Eq. (3.47 1 , it is clear that Eq. ( 3.50 1 may be taken to include r = if we replace VJ}^ /rhy t^l'^g-^'^ . 
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In particular, for a scalar field: 

(D' + r+ 1) Vr + l-Vr (e'„, - 4) + -^(D' - - iR')Vr-i = 0. (3.51) 



Together with the earlier equations, the transport equation Eq. (3.461 allows us to immediately solve for V^'^^ 
along a geodesic. To obtain the higher order Vy we also need to determine □'T^'if . At first sight this appears to 
require integrating along a family of neighbouring geodesies but, in fact, again we can write transport equations for 
it. First we note the identity 

CT a^a'a'^ d ...a'^ ' ' ' ^ R a'^^a' a'T a'^...c' (3.52) 

where U^Bcd = dcA^Bd-d.A'^Bd + A^cdA'^Bc-A^CdA'^Bc- Working with y/^' = A-VSv^ab'^ differentiating 
Eq. (3.501 we obtain a transport equation for the first derivative of 



p' + r + l)(F/^';,0+^'a'V;■^^^o' + a"'7^^'c'a'o'K^^' + ^ (a- (Ai/2v;lf)) ^0. (3.53) 

As noted above, this equation also includes r = if we replace /r in this case by g'^'~^ : 

(i?' + l)(yO^^^,0+r'a'^>0^^^«'+C7"'7^^'c'a'a'V^,^^' + ^ (A-I/^P^ (aVV^')) ^, =0. (3.54) 

Repeating the process 



{D'+r+ l)(V;r^^ .yy) + r b'V,''" ;a'a' + T a'V,''" 

i-iAB' _q' 7-.i3' i-rAC 

;/3' 



I ^q't^S' T/AC I ^Ol'tpB' i/AC I c"' i/AB' n' 75/3' T>AC' 

+ (T It C'b'a'Vr ;a' + CT /C C'a'a' V,- ;6' + 4 Q':h' ;a' - CT /t'^ a'6'a' 



+ ^'^>'7^^'c'a'o'V;^^' +a"'7^^'c'a'„';.'K^^' + ^ (A-V2pi^'^, (AV^y/f )) ^^^^ ^ 0, (3.55) 

with the Vff^^ -^a'V equation given by the same replacement as above. 

Clearly this process may be repeated as many times as necessary. At each stage we require two more derivatives on 
V;!^' than on l//*^' but this may be obtained by a bootstrap process grounded by the equation which involves 

only the fundamental objects A^/^ and g^'-^ which we have explored above. As with our previous equations, while 
this process quickly becomes too tedious to follow by hand it is straightforward to programme. 

For example, to determine Vi for a scalar field we first need to solve the two transport equations 

{D' + l)(l/o;a') + r'a'Vb;c' + \ {^'^'^ (□' " ".^ - ^R') A^/^) = 0, (3.56) 

and 

{D' + l){Vo;a'b') + r'6'Vb;a'c' + a'Vo;c'V + 

+ r'a';fc'Vbx' - a^' R"' a'b'c'Vo-.d' + +\ i^'^'^ (□' " - ^B!) A^'A = 0. (3.57) 

In the next two sections we show how the above system of transport equations can be solved either as a series 
expansion or numerically. For sufficiently simple spacetimes, it is also be possible to find closed form solutions which 
provide a useful check on our results. 

IV. SEMI-RECURSIVE APPROACH TO COVARIANT EXPANSIONS 



In this section, we will investigate solutions to the transport equations of Sec. Ill in the form of covariant series 
expansions. The goal is to find covariant series expressions for the HaMiDeW coefficients. Several methods have been 
previously applied for doing such calculations, both by hand and using computer algebra [2^ [Ul [5^1 ISZl HSl |M1 1301 
EllMllMllMlinSllMlETllMllMllinjIlIlllS- However, this effort has been focused primarily on the calculation of 
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the diagonal coefficients. To the our knowledge only the work of Decanini and Folacci [TH], upon which our method 
is based, has been concerned with the off-diagonal coefficients. 

Before proceeding further, it is helpful to see how covariant expansions behave under certain operations. First, 
applying the operator D' to the covariant expansion of any bi-tensor Aa-^...a^a{---a' about the point x, we get 

D'Aa^...a„^a{-a'^ {x , x') = ^ —k aa^...a^a[-a'„ /3i---/3fc {x)a'^^ ■ ■ ■ Cr^^^/Cr"' 

fc=0 

^ / 

= E -J-f" ^a,-a^<-< (4.1) 



fc=0 



where the second line is obtained by applying Eq. (3.5 1 to the first line. In other words, applying D' to the fc-th term 
in the series is equivalent to multiplying that term by k. 

Next we consider applying the operator D to the covariant expansion of any bi-tensor ^ai - a™a'j - a' about the point 
X. In this case, there will also be a term involving the derivative of the series coefficient, giving 



DAa,...a^a[---a',Sx,x') = ^ — fc aa,...a^a[-a'^ /3i • • • ft (2^) C^'^' ' •• ^'^'aCr" + Oai -a^ai •••«'„ I3,-I3k;a{x)cr'^' ' ■ ■ Cr'^" (T 



k=Q 



(4.2) 



n 

k=Q 



We can also consider multiplication of covariant expansions. For any two tensors, and _B°^, with product 
= A°'^B'^ij, say, we can relate their covariant expansions by 



00 



n=0 



01... 



\n=0 / \n=0 



= E ^ E (t)^"^ ft-..^". ,...-,y' (4-3) 



n=0 k=0 



Finally, many of the equations derived in the previous section contain terms involving the Riemann tensor at a;', 
i?" y^i^i- As all other quantities are expanded about x rather than x', we will also need to rewrite these Riemann 
terms in terms of their expansion about x: 



^1 (a|b|/3;7l---7t) 
fc=0 

E|^^M.)> (4.4) 



where we follow Avramidi [16j in introducing the definition 



b (n) = (ai|6|Q2;Q3-a„)'^ ' ' ' 

^£%(„)a"^-..a"". (4.5) 



These four considerations will now allow us to rewrite the transport equations of Sec. Ill as recursion relations for 
the coefficients of the covariant expansions of the tensors involved. 
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A. Recursion relation for coefficients of the covariant expansion of 7° 



Rewriting Eq. (3.15) in terms of covariant expansions, we get 



fc=0 



+ (£ (^^% (.)) (E • • • = 0. (4.6) 

From this, the n-th term in the covariant series expansion of 7°^, can be written recursively in terms of products of 
lower order terms in the series with /C: 



a' ira' a' r\ 

7 6 (0) = b 7 ^ (1) = 
7\ („) = - (^) E I („-.)7^ (4.7) 

Many of the following recursion relations will make use of these coefficients. 

B. Recursion relation for coefficients of the covariant expansion of rj°'^, 
Since 7" ^ is the inverse of yy'^jj, , we have 

r'o.v\'^S'^\,. (4.8) 

Substituting in covariant expansion expressions for 7" ^ and r/"^, , we find that the n-th term in the covariant series 
expansion of 77*^^, is 

V\ (0) = V\ (1) = v\, („) ^ (fc)ff"c.'7"; („_,)• (4.9) 

fc=2 ^ ^ 

C. Recursion relation for coefficients of the covariant expansion of ^"j,, 



Writing Eq. (3.141 in terms of covariant series, it is immediately apparent that the n-th term in the covariant 



expansion of , , is 



^ 6' (0) — " b' ? 6' (1) ~ " 

n-2 / . 

n 



k=2 



(„) = " 9'\rb- (n) - E U r= 7"'„ („_.)^^' (,). (4.10) 



D. Recursion relation for coefficients of the covariant expansion of A";, 



Using Eq. (3.171, we can write an equation for the n-th order coefficient of the covariant expansion of A°{,. However, 
the expression involves the operator D acting on the covariant series expansion of 77°^, , so we will first need to find an 
expression for that. As discussed in the beginning of this section, the derivative in D will affect both the coefficient 
and the cr^'s. When acting on the cr°'s, it has the effect of multiplying the term by n as was previously the case with 
D' . When acting on the coefficient, it will add a derivative to it and increase the order of the term (since we will 
then be adding a fx"). We now appeal to the fact that the terms in the expansion of 77"^, consists solely of products of 
/C"jj This means that applying the following general rules when D acts on /C^j^ ^^j, we will get the desired result: 

" b (fe) — '^'^ b {k) ^ b (fc+i) 
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• When encountering compound expressions (i.e. consisting of more than a single /C"^ ^j.^), use the normal rules 
for differentiation (product rule, distributivity, etc.) 

Taking this into consideration and letting £)"*" signify the contribution at one higher order and signify the contri- 
bution that keeps the order the same, we can write the general n-th term in the covariant series expansion of -D?7°{,/ 
as 

(i?ryV)(„) = D+7j\, („_i) - D^'tj^, („) (4.11) 
It is then straightforward to write an expression for the n-th term in the covariant series expansion of A^^: 

n-2 . V 
k=0 ^ ^ 

E. Recursion relation for coefficients of the covariant expansion of Aabc 



We can rewrite Eq. (3.241 as 

{D' + l){Aabc.l\) + Rabo^p^j'^l^, = 0, (4.13) 



which when rewritten in terms of covariant series becomes 

n-2 



1 fn\ fn\ 

Aabc (fe) = - — tE [A^'^aba {k)T c (n-k) + E [ k) '^"'"^ '^^^^'^ " 
I — n V / I — n V / 



(n-k) 



(4.14) 



k=0 ^ ' k=0 

where we follow Avramidi [161 117] in defining 

T^abc (n) = Rab{ai\c\;a2---a„)'^"^ ' ' ' O""" (4-15) 



Alternatively, writing Eq. (3.241 directly in terms of covariant series, we get 

n-2 



Aabc (k) 

1 1, ~r~ 1 / 1, ~r~ 1 

fc=0 

which has the benefit of requiring half as much computation as the previous expression. 

F. Recursion relation for coefficients of the covariant expansion of Babe 



By Eq. (3.26), we can immediately write an equation for the coefficients of the covariant expansion of Ba 



be- 



B, 



abc (n) 



1 " / \ 

-T.{l)^-baik)V"c(n-k) (4.17) 



n ^ V 

fc=o 



G. Covariant expansion of ^ 



From Eq. ( 3.34[ ) we immediately obtain expressions for the coefficients of the covariant series of 

C(o)=0 C(i)=0 ^(„^^_i_^p'^,^^^ (4.18) 

H. Recursion relation for A^'^^ 

Since C = InA^/^, we can write 

A^/^D'C = £i'Ai/2. (4.19) 
This allows us to write down a recursive equation for the coefficients of the covariant series expansion of A^/^, 



k=1 



^E(:)^C(.)A;/^.., (4.20) 
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I. Recursion relation for A 

Given the covariant series expansion for A^/^, it is straightforward to calculate the covariant expansion of A^^/^, 
which is simply 



-^\n) "Odd, 



^-1/2 ^ 

1 AJ/? n even. 
^. (") 



(4.21) 



J. Covariant expansion of r 



Eq. (3.371 may be immediately written as a covariant series equation, 



(4.22) 



K. Covariant expansion of r' 



Eq. (3.381 may be immediately written as a covariant series equation. 



T{n) 



(4.23) 



L. Covariant expansion of covariant derivative at x' of a bi-scalar 

Let T{x,x') be a general bi-scalar. Writing T as a covariant series, 

oo oo 

and applying a covariant derivative at x', we get 

oo 

T-b' = E ^(n);''' 

oo 

J(Qi - a„)'^ 1 • • - a " V ^, 

oo 

oo n y \ 

EEu)(^(^+i))(-i) 

n=Ofc=0 ^ ^ 



l)V{n-k) 



where we have introduced the notation 



{T(„))(-k) = T(ai. 



(4.24) 



•"(n-fc)1(n-fc + l)' 



(4.25) 



(4.26) 



M. Covariant expansion of d'Alembertian at x' of a bi-scalar 



Let T(x, x') be a general bi-scalar as in the previous section. Applying (4.251 twice and taking care to include the 
term involving , we can then write the d'Alembertian, □'T(a;, x') at x' in terms of covariant series. 



fc=0 



= - E ( /,)(^;a(/c+i))(-i)??(«-fc) - E ( L 



k=l 



k)- 



(4.27) 
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N. Covariant expansion of VaA^^^ 



Applying Eq. (4.251 to the case T = we get 



n=0 fc=0 



fe)- 



(4.28) 



O. Covariant expansion of □'A^'^^ 



Applying Eq. (4.271 to the case T = we get 



k=Q 



fc=l 



1/2 

;a(n — /c) ' 



where ^(n) is the n-th term in the covariant series of the tensor defined in (3.241 



(4.29) 



P. Covariant expansion of Vo 



The transport equation for Vq, Eq. (3.461, can be written in the alternative form 



{D' + 1) l/o - Vot' + i(n' - - ei?)Ai/2 = 0. 
This equation is then easily written in terms of covariant expansion coefficients, 



(n) 



h (E (fc)^o(.)r('„_,) - I ((□'A^/^)(„) - (m^ + Ci?)A 



1/2 

(n) 



(4.30) 



(4.31) 



Q. Covariant expansion of Vr 



The transport equation for Vr, Eq. (3.49) can also be written in the alternative form 



(D' + r+ 1) Vr - VrT' +—(□'- TO^ - ^R)Vr-l = 0. 

2r 



Again, this is easily written in terms of covariant expansion coefficients, 

Vr (r; 



(") - + (E (t)^r (fe)T('„_fe) - ^ ((□y.-Ocn) - (m^ + ei?)K-l („)) 



(4.32) 



(4.33) 



Tables [I| and illustrate the performance of our implementation on a desktop computer (2.4GHz Quad Core Processor 
with 8GB RAM). 



V. NUMERICAL SOLUTION OF TRANSPORT EQUATIONS 

In this section, we describe the implementation of the numerical solution of the transport equations of Sec. |III[ We 
use the analytic results for a, A^/^, ga'' and Vq for a scalar field in Nariai spacetimc from Ref. |33| as a check on our 
numerical code. 

For the purposes of numerical calculations, the operator V acting on a general bi-tensor T" ■j^, can be written as 



V'T" 



(5.1) 
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xiaiviiuevv L^oerncieni, dn 


Calculation Time (seconds) 


Number of terms 


]V[em.ory Used (bytes) 


fflo 








16 


ai 


0.004 


2 


288 


0.2 





7 


2 984 


ffla 


0.016 


68 


36 976 


ai 


0.168 


787 


522 952 


as 


3.012 


10183 


7 993 792 


aa 


70.196 


141691 


128 298192 


a? 


1016.696 


2 069 538 


2123 985 816 



TABLE I: Calculation performance of our semi-recursive implementation of the Avramidi method for computing the coincident 
(diagonal) HaMiDeW coefficients. 



Order 


Calculation Time (seconds) 


Number of terms 


Memory Used (bytes) 








2 


528 


1 





2 


288 


2 





8 


3 056 


3 


0.004 


12 


4 800 


4 


0.007 


41 


18 400 


5 


0.016 


72 


34 328 


6 


0.024 


189 


96 920 


7 


0.048 


357 


193 120 


8 


0.084 


810 


464 240 


9 


0.132 


1568 


938 960 


10 


0.188 


3 290 


2 067 512 


11 


0.384 


6 350 


4164 256 


12 


0.692 


12 732 


8 689 208 


13 


1.308 


24 340 


17230 944 


14 


2.688 


47 291 


34 697 312 


15 


5.156 


89 397 


67881528 


16 


9.692 


169 900 


133 241688 


17 


19.985 


317417 


256127 816 


18 


39.582 


593 371 


494 810 408 


19 


76.724 


1096 634 


937 590 072 


20 


122.943 


2 023 297 


1 766 643 856 



TABLE II: Calculation performance of the Avramidi method for computing the covariant series expansion of Vo- 



where s is the affine parameter, F^,^ 
we make use of the fact that 



are the ChristofFel symbols at x' and is the four velocity at x' . Additionally 



(s' - s)u'' 



(5.2) 



which allows us to write Eqs. ([sTT}, pl3| , pl8| , pl9| , p!22| ), ^1^, p!29| , ^A5\ and 



(3.461 as a system of coupled, tensor ordinary differential equations. These equations all have the general form 



^J^'-b'... = ('T'A'^' + S'^'fcV. .. + s'C^';; - T\;- T-^,p,u''' -■■■ + T^'-^, ^^'^^' + ' ' ' , (5-3) 

where we have set s = without loss of generality and where A" = initially (i.e. at s' = 0). It is not necessarily 
true, however, that the derivative of A" is zero initially. This fact is important when considering initial data for 
the numerical scheme. 

Solving this system of equations along with the geodesic equations for the spacetime of interest will then yield a 
numerical value for Vq. Moreover, since V = Vq along a null geodesic, the transport equation for Vo will effectively 
give the full value of V on the light-cone. We have implemented this numerical integration scheme for geodesies in 
Nariai and Schwarzschild spacetimes using the Runge-Kutta-Fehlberg method in GSL [IT. The source code of our 
implementation is available online [20j . 
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A. Initial Conditions 



Numerical integration of the transport equations requires initial conditions for each of the bi-tensors involved. These 
initial conditions are readily obtained by considering the covariant series expansions of Vq, A^/^, j^,, 77"^, and g„'' 
and their covariant derivatives at x' . Initial conditions for all bi-tensors used for calculating Vq are given in Table 
where we list the transport equation for the bi-tensor, the bi-tensor itself and its initial value 
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Additionally, as is indicated in Eq. (5.3 1, many of the transport equations will contain terms involving (s') ^. These 



terms must obviously be treated with care in any numerical implementation. Without loss of generality, we set s = 0. 
Then, for the initial time step (s' — s), we require analytic expressions for 



lim {s')-^A"'-^;^ 



(5.4) 



which are then be used to numerically compute an accurate initial value for the derivative. This limit can be computed 
from the first order term in the covariant series of A" ^ , which is found most easily by considering the covariant 
series of its constituent bi-tens ors. For this reason, we list in Table III the limit as s' ^ of all required bi-tensors 



with their limit as s ^ 0. 



multiplied by (s')"^ In Table IV we Hst the terms (s')"^^"V 



for each transport equation involving (s') ^, along 



Equation Bitensor 



Initial Condition 



(s') ^ Initial Condition 



(3.111 



(3.131 



5\ 



-9 b 



(3.181 



2 pa ,,Q 
'3-" (q|6|c)" 



(3.191 



bac^ — 3-rt (a\b\c)^ 



(^3.20 



2 r/a 

"3-"- {c\b\d) 



1 rta a 2 r>a a 2 r/a a 

2-"- (c|b|ti;a)" 3-"- (a|b|d)ic" 3-"- (ci|6|c);ti" 



(3.211 



1 pa 1 pa 

"3-"- {c\b\d) 2 bed 



1 pa , 2 pa Q 

2-" {c|6|d;ci)" + 3-"- (a\b\d);c^ 



(3.221 



9a' 



9a 



(3.241 



(3.291 



9a 



9a \c' d' 







2 acd 



2 aac 



2 pfc 

'3-"- ac(d;a) ' 



(3.341 



AV2 



Vo 



TABLE III: Initial conditions for tensors used in the numerical calculation of Vo- 



B. Results 



The accuracy of our numerical code may be verified by comparing with the results of Ref. |43j , which gives analytic 
expressions for all of the bi-tensors used in this paper. In Figs. [T] and [2] we compare analytic and numerical expressions 
for A^/^ and Vq, respectively. We consider the null geodesic which starts at p = 0.5 and moves inwards to p = 0.25 
before turning around and going out to p = 0.789, where it reaches a caustic (i.e. the edge of the normal neighborhood). 
The affine parameter, s, has been scaled so that it is equal to the angle coordinate, 0. We find that the numerical 
results faithfully match the analytic solution up the the boundary of the normal neighborhood. For the case of A^/^, 
Fig. [l] the error remains less than one part in 10^® to within a short distance of the normal neighborhood boundary. 
The results for Vo{x,x') are less accurate, but nonetheless the relative error remains less 1%. 

In Fig. [3] we give plot calculated from our numerical code which indicates how A'^^'^ varies over the whole light-cone 
in Schwarzschild. We find that it remains close to its initial value of 1 far away from the caustic. As geodesies get close 
to the caustic, A^/^ grows and is eventually singular at the caustic. This is exactly as expected; A^/^ is a measure 
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Equation 


Terms involving fs') ^ 


Initial condition for (s') 


^ terms 


1 J.ii 1 







( J.iJ 1 







(3.181 




II II 1 1 \ 

~ C J,/ / ^ t; j.i(J 1 1 — t iCr /,/ 
^ a' b'c' ^ d' a'c' ^ c' a.' o' J 


2 pa Q 


(J. it) 1 






-n c''^ b'c'-^ b'^'c'c'-^ c''^ c'b') 


1 oa OL 1 na 




(|3.20| 


— (J ft — 
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TABLE IV: Initial conditions for transport equations required for the numerical calculation of Vq. 
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FIG. 1: Comparison of numerical and exact analytic calculations of as a function of the angle, (/>, along an orbiting null 
geodesic in Nariai spacetime. Left: The numerical calculation (blue dots) is a close match with the analytic expression (red 
line). Right: With parameters so that the code completes in 1 minute, the relative error is within 0.0001% up to the boundary 
of the normal neighborhood (at = tt) . 

of the strength of focusing of geodesies, where values greater than 1 correspond to focusing and values less than 1 
correspond to de-focusing. At the caustic, geodesies are focused to a point one would expect A^/^ to be singular 
there. 

In Fig. |4] we give a similar plot (again calculated from our numerical code) wich indicates how V{x,x') varies over 
the light-cone in Schwarzschild. In this case there is considerably more structure than was previously the case with 
A^/^. There is the expected singularity at the caustic. However, travelling along a geodesic, V{x,x') also becomes 
negative for a period before turning positive and eventually becoming singular at the caustic. 

The transport equations may also be applied to calculate Vr{x,x') along a timelike geodesic. In Fig. [s] we apply 
our numerical code to the calculation of Vq{x, x') along the timelike circular orbit at r = lOAf in Schwarzschild. 
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FIG. 2: Comparison of numerical and exact analytic calculations of Vo for a massless, minimally coupled scalar field as a 
function of the angle, <j), along an orbiting null geodesic m Nariai spacetime. Left: The numerical calculation (blue dots) is a 
close match with the analytic expression (red line). The coincidence value is V{x,x) = |(C — |)/? = — |, as expected. Right: 
With parameters so that the code completes in 1 minute, the relative error in the numerical calculation is within 1% up to the 
boundary of the normal neighborhood (at = tt) , 





I 



FIG. 3: (Colour online) A^''^ along the light-cone in Schwarzschild Spacetime. The point x at the vertex of the cone is fixed at 
r = lOM. A}^'^ increases along a geodesic up to the caustic where it is singular. 



VI. DISCUSSION 



Several of the covariant expansion expressions computed by our code using the Avramidi method have been pre- 
viously given in Ref. [TH], albeit to considerably lower order (for example, Ref. [TH] gives V{x,x') to order (cr°)^ 
compared to order (cr") here). Comparison between the two results gives exact agreement, providing a reassuring 
confirmation of the accuracy of both our expressions and those of Ref. ^]. Furthermore, several of the expansions 
not given by Decanini and Folacci may be compared with those found by Christensen fl2; 13^. Again, we have found 
that our code is in exact agreement with Christensen's results. 

In Sec. |V] we discussed a numerical implementation of the transport equation approach to the calculation of Vq- 
This implementation is capable of computing Vq for any spacetime, although we have chosen Nariai and Schwarzschild 
spacetimes as examples. The choice of Nariai spacetime has the benefit that an expression for Vq is known exactly 
[43j . This makes it possible to compare our numerical results with the analytic expressions to determine both the 
validity of the approach and the accuracy of the numerical calculation. Given parameters allowing the code to run in 
under a minute, we find that the numerical implementation is accurate to less than 1% out as far as the location of 
the singularity of Vq, at the caustic point (i.e. everywhere within the normal neighborhood). 



Our Mathematica implementation of the semi-recursive approach to covariant expansions (Sec. IV I is given as a 
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FIG. 4: (Colour online) V{x,x') for a massless, scalar field along the light-cone in Schwarzschild Spacetime. The point x (the 
vertex of the cone) is fixed at r = IQM. V{x, x') is initially, then, travelling along a geodesic, it goes negative for a period 
before turning positive and eventually becoming singular at the caustic. 




FIG. 5: (Colour online) Vo{x,x') (solid blue line) and A^''^ (dashed purple line) for a massless, scalar field along the timelike 
circular orbit at r = lOM in Schwarzschild Spacetime as a function of the angle, (j) through which the geodesic has passed. In 
the logarithmic plot, the absolute value of Vo{x,x') is plotted, since Vo{x,x') is negative between (f) ~ 1.7 and (jj ~ 2.2. The 
apparent discontinuities in the plot are therefore simply a result of Vo{x, x') passing through at these points. 

proof of principle. We believe than considerable improvement could be achieved, particularly in the limiting area 
of memory requirements. The expressions for the HaMiDeW coefficients as computed by our code are very general. 
However, they are not necessarily given as a minimal set. For example, the DeWitt coefficient 03 may be written as 
a sum of four terms, yet our code produces a sum of seven equivalent terms. It is quite possible, however, that a set 
of transformation rules could be produced to reduce our expression to a canonical form such as that of Ref. [15 . In 
fact, some progress in that direction has been made already with the Invar ^45, 46, software package, which is able to 
quickly canonicalize scalar invariants. As our code is already written in Mathematica [37] and has the ability to output 
into the xTensor [55] notation used by Invar [JSJ HB], an extension of the package to allow for the canonicalization of 
tensor invariants would allow our expressions to be immediately canonicalized with no further effort. 

In integrating the transport equations given in Sec. |III| along a specific geodesic, we are not limited to the normal 
neighborhood. The only difficulty arises at caustics, where some bi-tensors such as A^^"^ and Vq become singular. 
However, this is not an insurmountable problem. The singular components may be separated out and methods of 
complex analysis employed to integrate through the caustics, beyond which the bi-tensors once more become regular 
(but not necessarily real- valued) [2 . This is highlighted in Fig. pi where our plot of A^^^ and Vq extends outside the 
normal neighborhood, the boundary of which is at (/i « 1.25, where the first null geodesic re-intersects the orbit. It 
does not necessarily follow, however, that the Green function outside the normal neighborhood is given by this value 
for V{x,x'). Instead, one might expect to obtain the Green function by considering the sum of the contributions 
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obtained by integrating along all geodesies connecting x and x' (there will be a discrete number of such geodesies 
except at caustics). 
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